Resummation of QED Perturbation Series by Sequence Transformations 
and the Prediction of Perturbative Coefficients 
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We propose a method for the resummation of divergent perturbative expansions in quantum elec- 
trodynamics and related field theories. The method is based on a nonlinear sequence transformation 
and uses as input data only the numerical values of a finite number of perturbative coefficients. The 
results obtained in this way are for alternating series superior to those obtained using Pade approx- 
imants. The nonlinear sequence transformation fulfills an accuracy-through-order relation and can 
be used to predict perturbative coefficients. In many cases, these predictions are closer to available 
analytic results than predictions obtained using the Pade method. 
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I. INTRODUCTION 

Perturbation theory leads to the expansion of a phys- 
ical quantity V{g) in powers of the couphng g, 



Tl = 



9 



(1) 



The natural question arises as to how the power series 
on the right-hand side is related to the (necessarily finite) 
quantity on the left. It was pointed out in |^ that pertur- 
bation theory is unlikely to converge in any Lagrangian 
field theory. Generically, the asymptotic behavior of the 
perturbative coefficients is assumed to be of the form B] 



(2) 



where K , 7 and S are constants. S is related to the first 
coefiicient of the (3 function of the underlying theory. 

In view of the probable divergence of perturbation 
expansions in higher order, a number of prescriptions 
have been proposed both for the resummation of diver- 
gent perturbation series and for the prediction of higher- 
order perturbative coefficients. A very important method 
is the Borel summation procedure whose application to 
QED perturbation series is discussed in The Borel 

method, while being useful for the resummation of diver- 
gent series, cannot be used for the prediction of higher- 
order perturbative coefficients in an obvious way. 

In recent years. Fade approximants have become the 
standard tool to overcome problems with slowly conver- 
gent and divergent power series Fade approximants 
have also been used for the prediction of unknown per- 
turbative coefficients in quantum field theory 1^-^. The 
[l/m] Fade approximant to the quantity Vig) represented 
by the power series (|^) is the ratio of two polynomials 
Pi{g) and Qm{g) of degree I and m, respectively. 



[l/m]v{g) = 



Qm{g) 



po+pig 



pi 9 



The polynomials Pi{g) and Qm{g) are constructed so that 
the Taylor expansion of the Fade approximation agrees 
with the original input series Eq. up to terms of order 
I + m in g, 



V{g) [l/m]r{g) = 0(5'+™+^) 



(3) 



For the recursive computation of Fade approximants we 
use Wynn's epsilon algorithm |^ , which in the case of the 
power series produces Fade approximants according 
to £2^'' = [n + k/k]-p{g). Further details can be found in 
Ch. 4 of fy. 

In this Letter, we advocate a different resummation 
scheme. For an infinite series whose partial sums are s„ = 
X]J=o ^J' nonlinear (Weniger) sequence transforma- 
tion with initial element sq is defined as [see Eq. (8.4-4) 
of 0]: 



V (-1V (^^ (/3+i)n-l Sj 



n\ (/3 + j)r 



1 



(4) 



jj (/3 + n)„-i ttj+i 



where {a)m — T{a + m)/T{a) is a Fochhammer symbol. 
The shift parameter (3 is usually chosen as /3 = 1 , and this 
choice will be exclusively used here (see also [0). The 
power of the 5 transformation and related transforma- 
tions [e.g., the Levin transformation, Eq. (7.3-9) of [Q] is 
due to the fact that explicit estimates for the truncation 
error of the series are incorporated into the convergence 
acceleration or resummation process (see Ch. 8 of [pO[). 
Note that the 5 transformation (|^) has lead to numer- 
ically stable and remarkably accurate results 1 11 , 1^ in 



1 
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the resummation of the perturbative series of the quartic, 
sextic and octic anharmonic oscillator whose coefficients 
display a similar factorial pattern of divergence as the 
quantum field theoretic coefficients indicated in Eq. (0). 
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TABLE I. Resummation of the perturbation series (g) for ~ 10. Results are given in terms of the dimensionless function 
Sb = 10^ [{8Tv^)/{-e^B^gB)] Sb- Apparent convergence is indicated by underlining. 



[(n + l)/2l/K2l 



1 


10.476 


10.476 190 476 


2 


-243.492 


- 1.617 535 903 


3 


10 530.918 


4.627 654 271 


4 


-774 888.106 


- 1.401 288 801 


5 


8.674 647 x lO'' 


2.773 159 300 



(0) 



'1,so(3b)) 



d(l,^o(5B)) 



-2.222 222 222 
-1.617 535 903 
^.820 833 551 
^.588 575 814 
-0.864 617 071 



-2.222 222 222 
-1.617 535 903 
4).820 833 551 
4).659 817 926 
-0.733 843 307 



60 
61 
62 
63 
64 
65 
66 
67 



-3.652 

5.553 
-8.721 

1.414 
-2.365 

4.082 
-7.261 

1.330 



544 X 
434 X 
566 X 
066 X 
759 X 
125 X 
275 X 
921 X 



^0209 
^0214 
^0218 

10 
10 

10' 



222 
226 
231 



-0.920 
^0.400 
-0.918 
-0.411 
-0.915 
-0.421 
-0.913 
-0.430 



487 125 
319 939 
054 104 
140 364 
746 814 
331 007 
555 178 
946 630 



5.992 
1.385 
-4.131 
-8.500 
2.890 
5.272 
-2.050 
-3.296 



187 X 
114 X 
495 X 
694 X 
004 X 
267 X 
491 X 
170 X 



13 



10^ 
10^ 
10^ 
10 
10" 
10" 
10^5 

10^5 



-0.805 633 9 81 
-0.805 633 9 80 
-0.805 633 97 9 
-0.805 633 97 8 
-0.805 633 97 7 
-0.805 633 97 6 
-0.805 633 975 
-0.805 633 975 



exact 



-0.805 633 975 



-0.805 633 975 



-0.805 633 975 



-0.805 633 975 



We consider as a model problem the QED effective ac- 
tion in the presence of a constant background magnetic 
field for which the exact nonperturbative result can be 
expressed as a proper-time integral: 



ds f , 1 s I / m* , , , 
— <icoths----;>cxp(- — .5) . (5) 



Here, B is the magnetic field strength, and e is the el- 
ementary charge. The general result for arbitrary E 
and B field can be found in Eq. (3.49) in |13| and 
in Eq. (4-123) in [Q. The nonperturbative result for 
can be expanded in powers of the effective coupling 
(7b = e^B^ /m'^, which results in the divergent asymptotic 
series 



2e'B- 



OO 

■ 5B X! 5b : 

n=0 



The expansion coefficients 

(_!)"+! 4" 



(2n-h4)(2n-|-3)(2n-t-2) 



(6) 



(7) 



where B2n+i is a Bernoulli number, display an alternat- 
ing sign pattern and grow factorially in absolute magni- 
tude, 



the partial sums Snign) of the perturbation series (g) 
produces convergent results even for a coupling constant 
as large as gs = 10. In the third column of Table |, we 
display the sequence 



[0/0], [1/0], [1/1] 



of Pade approximants, which were computed using 
Wynn's epsilon algorithm j|]. With the help of the no- 
tation \x\ for the integral part of x, the elements of this 
sequence of Pade approximants can be written compactly 
as [[[(n -I- l)/2]]/[Jn/2]]] . Obviously, Pade approximants 
converge too slowly to the exact result to be numerically 
useful. The Levin d transformation defined in Eq. (7.3-9) 
in [|lO|, which is included because it is closely related to 
the 5 transformation (^) , fails to accomplish a resumma- 
tion of the perturbation series, as shown in the fourth 
column of Table |. 

So far, predictions for unknown perturbative coeffi- 
cients were usually obtained using Pade approximants. 
The accuracy-through-order relation (||) implies that the 
Taylor expansion of a Pade approximant reproduces all 
terms used for its construction. The next coefficient ob- 
tained in this way is usually interpreted as the prediction 
for the first unknown series coefficient (see, e.g., [|| ||). 
The 5 transformation (|^), when applied to the partial 
sums Vn{g) of the power series (0), fulfills the accuracy- 
through-order relation ||ll|] : 



(-!)"+! r(2n-f2) 



.2ri+4 



(l + 0(2-(-+^))) (8) ng) - 5^\l.V,{g)) 



0{g 



.n+2\ 



g~.Q. (9) 



as n ^ OO. The series differs from "usual" perturbation 
series in quantum field theory by the distinctive property 
that all perturbation theory coefficients are known. 

The numerical results in the fifth column of Table ^ 
show that the application of the 6 transformation (0) to 



Upon re-expansion of the 5 transform a prediction for 
the next higher-order term in the perturbation series may 
therefore be obtained. 

In Table || we compare predictions for the coefficients 
Cn of the perturbation series (^) obtained by re-expanding 



2 



the Pade approximants [|n/2]/|(n — l)/2]]] and the 
transforms s!^}_2{l, so{gB)) , which were computed from 
the partial sums so(5b), silffB), ■ ■ • , Sn-i(5B)- For higher 
orders of perturbation theory in particular, the Weniger 
transformation yields clearly the best results, whereas 
for low orders the improvement over Pade predictions 
is only gradual. For example, let us assume that for 
a particular problem only three coefhcients cq, ci and 



C2 are available and C3 should be estimated by a ra- 
tional approximant. Because of the accidental equal- 
ity [l/l]p((7) = (5^"^ (1, 7^0(9)) , the predictions for C3 ob- 
tained using the Pade scheme and the S transformation, 
are equal. Differences between the Pade predictions and 
those obtained using the S transformation start to accu- 
mulate in higher order. 



TABLE II. Prediction of perturbative coefficients for tlie power series (g). Results are given for the scaled dimensionless 
power series = [(87r^)/(— e^B^pB)] Sb- First column: order of perturbation theory. Second column: exact coefficients. 
Third and fourth column: predictions obtained by re-expanding Fade approximants and Weniger transforms, respectively. 



exact 



K2l/[(n-l)/2l] Si%{l,so{gB)) 



3 +0.107 744 107 +0.050 793 650 +0.050 793 650 

4 -0.785 419 025 -0.457 096 214 -0.537 632 214 

14 -2.181 588 772 x 10^^ 



15 +2.055 682 756 x 10" +2.0 49 236 087 x 10"^ 

16 -2.199 481 257 x lO" -2.19 4 962 521 x 10^^ -2.199 48 091 x 10^** 

24 -1.711 360 421 x lO^'^ -1.711 272 235 x lO^'' -1.711 360 421 x 10^^ 

25 +4.421 625 118 x 10^^ +4.421 484 513 x 10^^ +4.421 625 118 x 10^" 

26 -1.234 699 825 x 10*^ -1.234 674 716 x 10^^ -1.234 699 825 x 10*^ 



We now turn to the case of the uniform background 
electric field, for which the effective action reads iflSi 



/"ds r , Is 





exp 



This result can be derived from (j^) by the replacements 
B ^ \ E and the inclusion of the converging factor. With 
the convention = E'^/m'^ the divergent asymptotic 
series 



2e^E^ 



n=0 



5e ^ 0, 



is obtained. The expansion coefhcients 



+4 



(2n + 4)(2n + 3)(2n + 2) 



(10) 



(11) 



display a nonalternating sign pattern, but are equal in 
magnitude to the magnetic field case [cf. Eq. (^)]. For 
physical values of gE, i.e., for gE > 0, there is a cut in 
the complex plane, and the nonvanishing imaginary part 
for Se gives the pair-production rate. As is well known, 
resummation procedures for (nonalternating) divergent 
series usually fail when the coupling g assumes values on 
the cut in the complex plane |jl^. The Borel method 
fails because of the poles on the integration contour in 
the Borel integral 0. The S transformation and Pade 



approximations fail for reasons discussed in [l^ and [p^ , 
respectively. 

We now come to an important observation which to the 
best of our knowledge has not yet been addressed in the 
literature: the prediction of perturbative coefficients by 
nonlinear sequence transformations may even work if the 
resummation of the divergent series fails, i.e. if the cou- 
pling g lies on the cut. A general divergent series whose 
coefficients are nonalternating in sign, evaluated for pos- 
itive coupling, corresponds to a series with alternating 
coefficients, evaluated for negative coupling. Alternat- 
ing series can be resummed with the S transformation in 
many cases, and predictions for higher-order coefficients 
should therefore be possible for both the alternating and 
the nonalternating case. For example, the perturbative 
coefficients in Eqs. (0) and ( |ll|) differ only in the sign 
pattern, not in their magnitude. As shown in Table III, 
rational approximants to the series (||) and (|l^) produce, 
after the re-expansion in the coupling, the same predic- 
tions up to the different sign pattern. 

We would like to stress here that the resummation 
procedure and the prediction scheme presented in this 
Letter also work for higher-order terms in the derivative 
expansion of the QED effective action The resum- 
mation also works for the partition function for the zero- 
dimensional theory which is discussed in (p. 464) 
and is used in as a paradigmatic example for the 
divergence of perturbative expansions in quantum field 
theory. Results will be presented in detail elsewhere p6| . 
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An interesting and more "realistic" application is given 
by the /3 function of the Higgs boson coupling in the stan- 
dard electroweak model In the MS renormalization 
scheme, five coefficients of this (3 function are known. 
Using the first four coefficients, the "prediction" for the 
fifth coefficient (which is known) may be obtained and 
compared to the analytic result. Using the transforma- 
tion a prediction of /34 w 4.404 x 10^ is obtained 
which is closer to the analytic result of (3^ ~ 4.913 x 10^ 
than the predictions obtained using the [2/1] and [1/2] 
Pade approximants (these yield /34 « 3.969 x 10^ and 
/34 « 4.188 x 10^, respectively). The prediction for 
the unknown coefficient /Js obtained using is /^s ~ 
—3.938 X 10^ as compared to (3^ 
[2/2] Pade approximant. 



(0) 

3.756 X 10^ from the 



TABLE III. Prediction of perturbative coefficients 
c'„ for the electric background field (p^. Results 
are given for the scaled dimensionless power series 
S'^=[{S^^)/{e'E^g^)\S^. 



exact 



4°22(i.^o(gE)) 



14 
15 
16 

24 
25 
26 



2.181 588 X 10^^ 
2.055 682 X lO" 
2.199 481 X 10^** 

1.711 360 X 10^^ 
4.421 625 X lO^'^ 
1.234 699 X 10"^ 



2.181 5 74 X 10^^ 
2.055 6 78 X lO" 
2.199 48 X 10^^ 

1.711 360 X 10^^ 
4.421 625 X 10^^ 
1.234 699 X 10"^ 



For the /3 function of the scalar (j)^ theory the situation 
is similar to the Higgs boson case. Five coefficients are 
known analytically [|l9[ . Again, the prediction for the 
fifth coefficient obtained using the transformation (Jj""* 
(1251.3) is closer to the analytic result of 1424.3 than 
the predictions from the [2/1] and [1/2] Pade approxi- 
mants which yield values of 1133.5 and 1187.5, respec- 
tively. For the unknown sixth coefficient, a prediction of 
— 1.70 X 10"* is obtained using (Jg"' whereas the [2/2] Pade 
approximant yields —1.63 x 10^. 

We have shown that the 5 transformation (^) can be 
used to accomplish a resummation of alternating diver- 
gent perturbation series whose coefficients diverge facto- 
rially. In many cases, the 5 transforms converge faster 
to the nonperturbative result than Pade approximants. 
The 6 transformation uses as input data only the numer- 
ical values of a finite number of perturbative coefficients. 
We stress here that the factorial divergence is expected of 
general perturbative expansions in quantum field theory 
[see Eq. (|])]. The Weniger i5 transformation can be used 



for the prediction of higher-order coefficients of alternat- 
ing and nonalternating factorially divergent perturbation 
series. Both in model problems and in more realistic ap- 
plications, the 5 transformation yields improved predic- 
tions (compared to Pade approximants). It appears that 
the potential of sequence transformations, notably the 5 
transformation, has not yet been widely noticed in the 
field of large-order perturbation theory. 
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